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ABSTRACT 

The stochastic quantization of dissipative systems is discussed. It is shown that 
in order to stochastically quantize a system with dissipation, one has to restrict 
the Fourier transform of the space-time variable to the positive half domain in the 
complex plane. This breaks the time-reversal invariance, which manifests in the 
formulation through the resulting noninvariant forms for the propagators. The 
relation of the stochastic approach with the Caldeira and Leggett path-integral 
method is also analyzed. 



f tzani@vana.physto.se; address after September 1994: Imperial College, Physics Dept., 
Prince Consort Road, London SW72BZ 



The quantization of a classical system may, in general, be hindered by sev- 
eral obstacles, such as nonrenormalizability, topological obstructions, anomalies 
etc. There exists, moreover, a whole class of systems for which there is no avail- 
able method of quantization, namely non-hamiltonian systems. The quantization 
of systems with no lagrangian or hamiltonian formulation, apart from the well- 
known examples of systems with dissipation, includes several interesting physical 
problems, such as the Bargmann-Wigner higher spin field equations [1] and field 
theories of fundamental anyons [2]. For the Bargmann-Wigner system of equa- 
tions there is no action which directly reproduces this set of equations, although 
the quantization for higher spin systems has been done using more complicated 
and indirect ways [3]. The case of anyonic systems is known to be unsolved. The 
most common example of a non-hamiltonian system is, nevertheless, the case of a 
particle on the line with friction. Its equation of motion is 

■ n 

mx + 7X H — ;— = U (1) 
ax 

where overdot denotes time derivative and V{x) is the potential of the particle. 
If we assume an action of the form 

S = J f{x,x,t)dt (2) 

for the above equation, the variation of the action gives 

^= fdt(^-^^-^^x-^x) (3) 
5x J dx dt dx dx dx dxP- 

Assuming, now, that the second derivative term of the equation is multiplied by 
a constant, that is, -q^ — m, the function / can be written as f — f\{x,t) -\- 
f2{x,t)x + ^miP'. Putting this form of / into (3) we obtain 



5S f dh df2 ... 



which never contains first derivative terms. Therefore, there is no action from 
the variation of which the friction term can be derived and therefore no canonical 
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structure and no quantum mechanics for the system (1). (Notice that in the search 
for an action above we do not allow for an equation of motion which is multiplied 
by a time- dependent factor. Such time-dependent factors of the form e-^*-*-* have 
been considered [4] in the past, but they lead to a time-dependent hamiltonian and 
therefore not desirable physics.) 

The model (1) (and its higher dimensional versions), apart from its interest as 
a phenomenological description of dissipation due to interactions, also has recently 
found applications in string theory [5]. 

There are mainly two approaches that have been used in the literature in the 
study of this problem. The first consists of considering the quantization of these 
systems as a fundamental physical question and try to find a direct quantization 
procedure which classically reduces to the known dynamics [6,7]. The second con- 
sists in coupling the system with a many degrees of freedom reservoir. The system 
plus the reservoir, now, can be described by a hamiltonian and the dissipation is 
due to the interaction of the system with the reservoir. In this case the problem 
is reduced to finding a way to "integrate out" the extra degrees of freedom of the 
reservoir in order to finally obtain a quantum mechanical description of the original 
system [8-9] . The first approach is more fundamental and mathematical while the 
second is more physical. 

Using the second approach in the study of tunneling in dissipative systems, 
Caldeira and Leggett [9] have derived an effective action for (1) in the euclidean 
time. In their path-integral approach the term which classically produces the 
friction force corresponds to a non-local term in the action. Their effective action 
is given by 

SieffMt)) = J ilmx' + V{x)]dt + ^ j j dtdA^^^^0^ (4) 

-oo 

The approach of Caldeira and Leggett has, however, some disadvantages. Coupling 
the system to an infinite set of oscillators and then integrating them out, leads to a 
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state where the fourth moment of the field variable diverges. This is due to the zero 
point motion of the set of oscillators, which perturb the system in a substantial way 
[10]. Moreover, it is not clear that this approach is equivalent to a fundamental 
quantization procedure of dissipative systems. Nevertheless, this action has been 
used in the literature in a similar context in order to describe the quantum effect 
of friction [11]. 

The first approach was followed by Kostin [6] who has proposed a modified 
Schrodinger equation for the model (1) in which the friction is reproduced through 
a wavefunction-dependent potential. This approach exhibits also some drowbacks: 
it violates the superposition principle, and has stationary states in which the energy 
does not dissipate. (The same equation was rederived by Yasue [12] using Nelson's 
stochastic quantization scheme [13].) Although we do not expect the quantum 
mechanics of dissipative systems to have the same properties as the standard ones, 
the fact that one obtains stationary non-dissipative states for these systems is 
counter-intutive. Due to these rather unphysical properties Kostin's equation has 
not been used extensively in the literature. 

This last approach of a modified Schrodinger equation has been extended fur- 
ther by Polychronakos and the author in a previous paper [10] by allowing a 
wavefunction-dependent first-derivative term. This study resulted in a paramet- 
ric non-linear Schrodinger equation for the particle with friction which reproduces 
Kostin's equation for a specific value of the parameter. This family of equations 
avoids some of the problems that Kostin's equation exhibited, such as the station- 
ary states. 

In the present paper, taking the point of view of fundamentally quantizing the 
friction system, we study the stochastic quantization approach of equation (1). 
Stochastic quantization [14] has the advantage that it relies only on the classi- 
cal equation of motion. It essentially consists of reproducing the weighting factor 
exp(— 5') of the euclidean path- integral as a limiting probability density of a partic- 
ular stochastic process. The stochastic process which reproduces the path-integral 
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results for bosonic fields is given by the following Langevin equation 

where ^ is the classical equation of motion, r is the stochastic time and is a 
white noise over spacetime and r. In the stochastic limit r — > oo the fields attain 
their equilibrium distribution with probability density equal to exp{—S). 

Since the stochastic quantization method is based on the equation of motion, 
one would naively expect to overcome the problem of hamiltonian formulation of 
these systems and be able to quantize them directly using stochastic quantization 
method. In this work, we presend some peculiarities related with the stochastic 
quantization approach of dissipative systems. 

We consider the simplest example of a particle moving on the line with friction. 

Its equation of motion is given by (1) 

mx + 272; + kx = 0, (6) 

where we have multiplied the friction term by 2 for later convenience and have 
assumed the harmonic oscillator potential V{x) = ^kx^. The Langevin corre- 
sponding to equation (6) is given by the following euclidean time equation 

dx{t, t) ^ _ 2 -^^^^^ _ ^^^^^ ^ ^^^^ (j^ 



dr 



where 



< ri{t, r) >^ = and < rj{t, T)rj{t', t') >^ = 25{t - t')5{T - r') . 

The first thing to notice about this Langevin is that it is complex. Since the 
friction term is first order in time derivative, after the Wick rotation to euclidean 
time this term is imaginary in the expression of Therefore, one faces the 
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problem of stochastically quantizing complex actions [15] which has not been solved 
in its generality. In this particular case, however, since the equation of motion is 
linear in the x variable, the problem can be solved. The solution of the Langevin 
equation is easily obtained in the Fourier space. The Fourier transform of (7) is 

dx{u, t) ^ ^_^2 ^ _ k)x{uj, t) + fi{uj, t) (8) 

UT 

where x{u},t) and fj{uj,T) are the Fourier transforms of x{t,T) and rj{t,T) corre- 
spondingly. (Notice that for overcritical damping, that is 7 > 7c = \/km, (8) will 
develop instabilities for values of a; such that ^ — \/^ — -^ < u < ^-\-\^ — ■^.) 

Now, in the Fourier space the Langevin becomes real and it naively seems that 
one avoids the complex action problem by going to Fourier space. This apparent 
resolution is illusionary, however, and is due to the Fourier transform of first order 
derivatives (?7(c<j, r) is complex). As we will see in what follows, the problem related 
with the first order in time derivative term manifests itself in the impossibility to 
obtain 7-dependent results in this approach. The solution of the last equation is 
given by 



where now the ?7-average in the Fourier space is 

< 77(0;, r) 77(0;!, ri) >^=25(a; + a;i)5(T-ri). (10) 

The interesting quantities to be computed are correlation functions of the vari- 
able x{i) . In the Fourier space they are 

< x{u),T)x{0Ji,Ti) >rj 



T ri 







6 



(/c — 27a; + mcj^j + (/e — 270;! + moj^) 

where we have used the relation (10) and have taken the hmit r = ri — > 00 in 
order to obtain the last result. 

Next, we Fourier transform the last two-point function back to the original 

space. Because of the form of the white noise average in the Fourier space, namely, 
given as a delta function of the sum of the c<;'s in (10), all linear terms in to cancel 
after the integration over lo. We obtain 

< x{t)x{t') >= 7^ / dio i-— e^'^(*-*') (12) 



This integral is calculated by contour integration. The poles {u = ±zy — ) lie on 
the imaginary axis and for t > t' we close the contour in the upper half complex 
plane while for t < t' we close the contour in the lower half plane. The result 
obtained is 

< x(t)x(t') > = -^=[9it - t')e-^^'-*'^ + e{t' - t)ev^(*-^')] (13) 
2\/km 



which gives 



< x\t) > = -= (14) 
y km 



for the equal time correlation function of the variable x{t). These results coincide 
with the results one obtains from the quantization of this system for the 7 = 
case. Therefore, the naive application of stochastic quantization for this system 
gives results for the correlation functions that do not depend on the friction term. 
(It is a straightforward calculation to show that none of the higher order correlation 
functions depend on the friction parameter.) 

The meaning of the last result is the following: adding a friction term in the 
equation of motion is equivalent to adding a "curl" term in the Fokker- Planck 
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hamiltonian. Indeed, the Fokker-Planck equation for the system is 



dt 



5 f 5 5Sq 



5x{t) \5x(t) 5x(t) 



dP 



where is the part of the equation of motion which is derived from an action 
and F is the friction force, given by —2i'yx{t) in our case. Because of the form 
of the Fokker-Planck equation we can identify the probabihty current as Jx — 
— + ^^t) +-^)^ is probabihty density, which satisfy the probabihty 

conservation equation 

6Jx dP 
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Then one can show that the part of the current Jf which comes from the friction 
force is divergenceless. Indeed, the Fokker-Planck equation in the equilibrium limit 
is satisfied by P = e~^°, and 

(The last relation is easily shown in the u! plane and it is based on the fact that 
the range of integration in the uj variable is from — oo to oo while the integrand is 
an odd function of oo.) That is, V ■ = 0. Therefore, the friction force produces 
currents that do not change the probability density, which corresponds to the fact 
that they do not alter the physical results. 

One, however, can approach the problem in the following way: let's view x{t) 
as a -|- 1-dimensional field theory. Then, for the case 7 = 0, x{t) is real and 
it describes only particles. (Particles here correspond to energy quanta.) In this 
description, in the Feynman picture, positive frequency uj > which propagates 
forward in time and negative frequency a; < which propagates backward in time 
corresponds to the same particle state. Since friction breaks the time-reversal 
invariance, this is not a good description for the 7 7^ case. There exists, however, 
an equivalent description for particles in the frictionless case. This can be obtained 
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by considering a complex field but keeping only positive frequencies and forward 
propagation in time. This is, then, a more appropriate way to incorporate friction. 
In order to achieve it, we make the assumption that the field x{t) in our Langevin 
equation is complex with only positive frequencies, that is, an analytic function in 
t in the upper half complex plane. 

Given a real function x{t) we can always construct a function analytic in the 
upper half plane (in general complex), with the use of the Hilbert transform of the 
function, as follows 



where by P in front of the integral we mean principal value. With this Xa{t), 
Xa{co) — for any a; < 0. We assume that the Langevin equation (7) is satisfied by 
the function Xa{t,T). Then, in order for the Langevin to be consistent, the white 
noise function should obey the same analyticity property. That is 



In order to calculate correlation functions in this framework we must remember 
that our Langevin describes only particles and not antiparticles. Then, for t > t' 
the appropriate correlation function to be computed is 



This gives the propagator for the particles moving forward in time from t' to t. 
(The propagator for the antiparticles is given by the complex conjugate of this 

quantity.) 

Again it is convenient to go to the Fourier space. The Langevin is given 
by (8) where x{u!,r) and fj{u!,r) have been replaced by Xa{uj,r) and 770(0;, r) 
correspondingly. Since the white noise is an analytic function in the t-plane it 




(18) 




(19) 



< Xa{t)x*^it') > 



(20) 
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satisfies the condition r]a{uj,T) = for any c<j < in tlie a;-piane. Tlien tlie 
correiation functions for fja obtained from (10) are 

<fia{u;,T)fia{uj',T') > = 

(21) 

< U^, r)fj*^{u;', r') > = 8 S(uj - u;')5(r - r') 

for botli LU and oj' greater than zero. 

Following the same procedure as before and using the last relations for the rj- 
average we obtain the following expression for the two-point function of the analytic 

Xa{t) 

oo 

<x.W<(0>-^/^^ „,._;^,^, eM.-0 (22) 



which depends explicitly on the friction parameter. 

Next, we compute the last integral. We are interested only in the case that 
t > t'. For t > t' wc can evaluate the integral by deforming the contour into the 
upper half uj plane. The integrand has poles at u = ^ ±ia where a = \J ^ 
and we have assumed that 7^ < km . If we further make the assumption that we 
are close to critical damping, that is, 7^ — km is a very small quantity, the poles will 
he close to the real axis. We deform the contour as shown in fig. 1. Apphcation 
of Cauchy's theorem around the contour C gives the following expression for the 
above integral 

e+ioo 

4 /■ , 1 



27r 



/ du^^ ^e-(*-*') + -^e(^-")(*-*') (23) 

J mijj''- — 27a; + k 2ma 



Then if we assume that \t\^ » 1 and \t\a < 1 we can show that the integral 
in the above expression is negligible and we obtain the following result 

< xa{t)xl{t') >= ^ e(^^-\/^-^)(^-^') (24) 
Y km — 7^ 

The last expression gives the propagator for the problem of a particle with 
friction. In Minkowski space this propagator is not unitary due to the reality of 



10 



the term proportional to ^ in the exponent. This means that there is particle 
loss in this picture, which corresponds to the decay of energy due to dissipation, 
since particles correspond to energy quanta. Indeed, from the classical equation of 
motion, the change in the energy (E) is given by 

E = -27i;2 (25) 

Choosing, now, our state to be an almost "energy eigenstate" (which can be the 
case for 7 << 1) we obtain 

< E >^ -2^ <x'^ >^ < E > (26) 

m 

for the decay of the energy due to the friction. This is the same as the decay rate 
we obtain for the energy from the result (24). 

Finally, one could solve the Langevin equation (7) by separating its real and 
imaginary parts. We obtain the following coupled equations for x{t) and y{t). 

^^^^ = mx{t, t) + 2^y{t, r) - kx{t, r) + r}{t, r) 

and 

dv(t t) 

= my{t, t) - 2-fx{t, r) - ky{t, r) + r]i{t, r) (27) 

where y{t) — —^P J jjz^dt' and y{t) — |P / jjzfjsdt' and by principal value, 
here, we mean 



J {t- t'Y 2 [7 {t-t'^ ieY J {t-f- iey 
The stochastic process is, then, defined by the equations (27) . 



(28) 
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The corresponding equation for x{u,t) is 

^^^^'^^ = i-mw'^ + 2^\uj\ - k)x{u;, r) + fj{uj, r) (29) 
or 

The last equation coincides with the original Langevin (8) after the substitution 
of a; by I a; I and gives results for the correlation functions of x{t) which depend on 
the friction parameter. 

Indeed, with this Langevin the two-point function is expressed as an integral 
over u) in the following way 

00 

<.W.(0>=^/d. ^,_;,^,^^ e-('-'-) (30) 

— OO 

The calculation of this integral follows similar lines as the previous one and 
the result for the correlation function is given by ^ of the real part of the expres- 
sion (24). (Remember that x{t) here is the real part of ^^(t) and that y{uj) = 
—isgn{uj)x{uj), where by sgn{uj) wc denote the sign of uo.) This expression for 
7 = reproduces the standard harmonic oscillator result (12). The full expression 
(24) can be recovered from (30) by using dispersion relations. 

In what follows, we compare the stochastic approach with the path-integral 
approach of Caldeira and Leggett. We find that the stochastic quantization one 
obtains starting from the Langevin given by equations (27) is equivalent to the 

stochastic quantization one would obtain starting from the effective action of 
Caldeira and Leggett. 

In order to see this, it is convenient to rewrite the non-local term of the action 
(4) with the use of the identity x{t)x{t') = ^[x'^(t) + x^(t')] - l[x(t) - x(t')f. The 
effective action can be expressed as 

r OO r oo r 

-oo -oo 

(31) 

where the principal value in the last integrals is necessary in order to take care 
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of the infinity at t = t' . The last modification one has to do is to extend the 
hmits of integration in (3f ) from — oo to oo [9] . It can be easily, then, shown that 
the last integral vanishes and the variation with respect to x{t) of the remaining 
7-dependent term in (31) gives 



which coincides with the friction term of the first of the equations (27), was not 
for the factor of 2 we have multiplied our friction term. 

Concluding, we like to emphasize that it seems peculiar that the stochastic 
quantization of our dissipative model is achieved by restricting the Fourier trans- 
form of the space-time variable to the positive half domain in the complex plane. 
The physics, however, dictates this choise as it is argued in the text. Choosing the 
solution of the Langevin equation to be a complex function is natural, since in our 
original Langevin x{t) was complex. The extra restriction, however, to be analytic 
in t in the upper half plane stems from the fact that the Langevin should break 
the time-reversal invariance in order to describe friction. Notice at this point that 
if one complexifies both x{t) and "qit) in the Langevin one obtains doubling in the 
degrees of freedom of the problem. This doubling is eliminated after restricting 
the frequencies in the positive half domain in the complex plane. 

It is interesting, nevertheless, to connect our approach with what is known 
about this problem. Due to the time-reversal non-invariance of the dissipative pro- 
cesses, retarded Green's functions are more appropriate to use in order to describe 
them than the symmetric time-ordered ones [16]. The reason is that, since the 
system looses energy irreversibly, a function with only forward propagation in time 
can describe it. The standard stochastic approach, however, gives the Feynman 
propagator, which is symmetric under time reversal. Indeed, our two-point func- 
tion (13) obtained by straightforward application of stochastic quantization does 
not break time-reversal and does not describe dissipation. Formulating, on the 
other hand, the problem such that it involves fields which are analytic in the upper 




(32) 
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half plane (or equivalently fields that contain only positive frequencies) we have 
broken time- reversal in such a way that when time runs from — oo to oo these fields 
are appropriate to describe the forward propagation in time. The complex conju- 
gate of these fields gives the antiparticle propagation, which can be equivalently 
obtained by performing a time reversal operation in the original fields. The exact 
relation, however, of our approach with the retarded Green's function approach, or 
equivalently with the Keldysh formulation [16-17] for nonequilibrium phenomena 
in quantum field theory, is a subject for further investigation. 

It is also interesting to notice that the fact that the stochastic approach with 
the uj restricted in the positive half plane gives friction-dependent results is related 
to the divergence of the corresponding probability current as explained earlier. 
Indeed, in this case, due to the fact that the range of integration of lo is from to 
oo, the calculation of the divergence of the J/ current leads to a nonzero result. 
Therefore, the friction currents change the probability density in this case, as they 
should. 

Finally, in general, the stochastic quantization of dissipative systems will work 
in similar patterns; that is, one has to complexify each real coordinate of the system 
and keep only positive frequencies. 

The material of this paper was presented at the "Lattice Field Theory Work- 
shop", Vienna, June 11-12, 1993 and at the "IV International conference on math- 
ematical physiscs, string theory and quantum gravity" , Alushta, 13-24 June, 1993. 

After the completion of this work, preprint [18] appeared where the prescription 
a; — > I a; I is used in the finite temperature stochastic quantization of the system. 
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